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ABSTRACT 

We present the analysis of a new broad-band X-ray spectrum (0.6 — 180.0 keV) of 
the radio-quiet quasar MR 2251-178 which uses data obtained with both Suzaku and 
the Swift /HAT. In accordance with previous observations, we find that the general 
continuum can be well described by a power-law with F = 1.6 and an apparent 
soft-excess below 1 keV. Warm absorption is clearly present and absorption lines due 
to the FeUTA, FeL (Fexxiii-xxiv), Sxv and Sxvi are detected below 3keV. At 
higher energies, FeK absorption from Fexxv-xxvi is detected and a relatively weak 
{EW = 25tl^ eV) narrow FeKa emission line is observed at E = 6.44 ± 0.04 keV. The 
Fe Ka emission is well modelled by the presence of a mildly ionised < 30) reflection 
component with a low reflection fraction (i? < 0.2). At least 5 ionised absorption 
components with 10^° < A'h ^ 10^"^ cm^^ and 0< logf/ergcms"^ <4 are required to 
achieve an adequate spectral fit. Alternatively, we show that the continuum can also 
be fit if a F 2.0 power-law is absorbed by a column of Nn ~ 10^'^ cm~^ which covers 
~ 30% of the source flux. Independent of which continuum model is adopted, the FeL 
and Fexxv Hea lines are described by a single absorber outflowing with Wout ~ 0.14 c. 
Such an outflow/disk wind is likely to be substantially clumped {b ^ 10^'^) in order 
to not vastly exceed the likely accretion rate of the source. 
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1 INTRODUCTION 

It is now well established that the soft X-ray spectrum of 
at least half of all Seyfert 1 galaxies is characterised by re- 
gions of photoionised "warm" absorption along the line of 
sight (Blustin et al. 2005; McKernan et al. 2007). Extensive 
studies of the local Seyfert population with Chandra and 
XMM-Newton have revealed that the warm absorber (WA) 
is an intricate array of narrow absorption lines, often seen 
to be moderately blueshifted with v ~ 100 — 1000 kms~^, 
from various ionisation stages of abundant elements such as 
C, N, O, Ne, Mg, Si, S and Fe (e.g. Kaastra et al. 2000; 
Kaspi et al. 2002; Crenshaw, Kraemer & George 2003). De- 



tailed modelling of the WA with photoionisation codes such 
as XSTAR (Bautista & Kallman 2001) have shown that the 
absorbing material is typified by an ionisation parameter in 
the range log^/ergcms~^ ~ — 3, a column density of be- 
tween A^H ~ lO'^" — 10^'' cm~^, and most likely originate in 
a wind outflowing from the putative torus (Blustin et al. 
2005) or the later stages of an accretion disk wind (Proga 
& Kallman 2004). Energetically, due to their low Vout, soft 
X-ray WAs typically have kinetic luminosities to the order 
of ~ 1% of an AGN's bolometric luminosity (e.g. Blustin et 
al. 2005) and are therefore unlikely play an important role 
in terms of AGN feedback scenarios. 

In addition, there is now also a large body of obser- 
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vational evidence for blueshifted absorption lines at rest- 
frame energies greater than 7keV in many AGN (e.g. 
PG 1211+143, Pounds et al. 2003; PDS456, Reeves et al. 
2003; NGC1365, Risaliti et al. 2005; MCG -5-23-16, Braito 
et al. 2007; H 1413-f 117, Chartas et al. 2007; Mrk766, Miller 
et al. 2007; NGC 3516, Turner et al. 2008; Mrk 509, Cappi 
et al. 2009; 30 445, Reeves et al. 2010 & Braito et al. 2011; 
Tombesi et al. 2010a; Tombesi et al. 2010b). This absorp- 
tion is generally attributed to Fexxv/xxvi in an outflowing 
accretion disk wind and requires both a high column den- 
sity (i.e. A'^H = 10^^ — lO^** cm~^ or above), high ionisation 
parameter (i.e. log^/ergcms~^ ~ 3 — 6), and an outflow 
velocity ~ 0.1 c. 

Given their large outflow velocities the resulting kinetic 
power of accretion disk winds can reach a significant fraction 
of a sources bolometric luminosity (e.g. Pounds & Reeves 
2007; Tombesi et al. 2010b) and they may represent a driving 
mechanism for AGN feedback processes. Indeed, extrapolat- 
ing the kinetic luminosity over a conservative active phase 
and outflow duty cycle for a given AGN, it is often found 
that the total mechanical output can be in excess of, or is 
at least very similar to, the binding energy of a reasonably 
sized galaxy bulge. Accretion disk winds could therefore be 
a possible explanation for the observed M — a relation (Fer- 
rarese & Merritt 2000; Gebhardt et al. 2000); where mass 
ejected from a radiatively driven accretion disk wind (e.g. 
King 2003) can move out into a host galaxy, sweeping up ma- 
terial until, once a critical black hole mass has been reached, 
the ISM has been evacuated of material and both star forma- 
tion and SMBH growth cease (see King 2010, and references 
therein) 

1.1 The Radio-Quiet Quasar MR 2251-178 

MR 2251-178 (z = 0.064; Bergeron et al. 1983; Canizares et 
al. 1978) was first detected as a bright X-ray source during 
the Ariel V all-sky survey (Cooke et al. 1978) and sub- 
sequently identified as a radio-quiet quasar by Ricker et 
al. (1978) using SAS-3 observations. The quasar is located 
on the outskirts of a cluster of approximately 50 galaxies 
(Phillips 1980) and is surrounded by a large extended neb- 
ula of diffuse gas which is characterised by [O ill] emission 
in the optical (Bergeron et al. 1983). More recent observa- 
tions have also noted similar extended emission in the X-ray 
band (Kaspi et al. 2004; Gibson et al. 2005). The source is 
observed to be a weak radio emitter, and has a Fanaroff- 
Riley type I (FRI) morphology. 

The first detailed spectral study of MR 2251-178 in the 
X-rays was conducted by Halpern (1984) who noticed that 
absorption in the soft X-ray band varied on time scales of 
around 1 year. The variability implied an order of magnitude 
increase in the absorption column and was attributed to 
the change in ionisation of material along the line of sight. 
Historically, this is regarded as the first suggestion of an 
X-ray AGN warm absorber. 

Subsequent observations with EXOSAT, Ginga and 
BeppoSAX established that the broadband X-ray spectrum 
of MR 2251-178 can be well described by a power-law of 
photon-index F ~ 1.7 which is absorbed by a column den- 
sity of around a few xlO^^ cm~^ (Pan et al. 1990; Mineo & 
Stewart 1993) , and a high-energy roll-over at around 100 keV 
(Orr et al. 2001). Mineo & Stewart (1993) also found that 



the ionisation state of the absorbing material was strongly 
correlated with the source luminosity. 

In the ultra-violet (UV), Monier et al. (2001) found ab- 
sorption lines due to Lya, Nv and Civ blueshifted in the 
rest-frame by ~ 300kms~^. The Civ absorption was later 
shown by Ganguly et al. (2001) to vary over a period of 
roughly 4 years which enabled the authors to deduce a max- 
imum distance of ~2.4kpc between the absorption clouds 
and the continuum source. 

Kaspi et al. (2004), using a series of ASCA, BeppoSAX 
and XMM-Newton observations which spanned 8.5 years, 
confirmed that the continuum can be described by an ab- 
sorbed power-law of photon index F ~ 1.6 but found that it 
further required an additional soft excess at low X-ray en- 
ergies. The XMM-Newton spectrum required at least two or 
three separate ionised absorbers with column densities in the 
range IQ-^""^-^ cm~^ and had physical properties which ap- 
peared to vary between observations. This lead the authors 
to posit a scenario where absorption clouds were moving 
across the line of sight over the timescale of 'several months'. 
Further UV absorption lines from C ill, H I and O vi were de- 
tected in the FUSE spectrum, which were blueshifted with 
velocities similar to those found by Monier et al. (2001). 

MR 2251-178 has also been observed by the Chan- 
dra/HETG which revealed evidence of a highly-ionised, 
high-velocity (i.e. v = -12700 ± 2400kms~^) outflow in 
the FeK band (Gibson et al. 2005). Gibson et al. (2005) 
attributed this absorption to the Lya line of Fexxvi, and 
inferred that unless the absorber has a low covering fraction 
the mass-loss rate of MR 2251-178 is at least an order of 
magnitude larger than the accretion rate. 

In this paper we present a new broad-band X-ray ob- 
servation of MR 2251-178, using data obtained with both 
Suzaku and the Swift Burst Alert Telescope (BAT). We 
begin by discussing the data reduction process in iJ2I a^nd 
then parameterise the broadband spectrum in 331 In >i3.2l 
we discuss the several absorption lines which are detected 
and determine their statistical significance, before perform- 
ing a detailed modelling of the spectrum with the xstar 
photoionisation code in ^ To calculate luminosities, a con- 
cordance cosmology with Ho = 71 kms~^ Mpc~^, f2A=0.73, 
and nm=0.27 (Spergel et al. 2003) was adopted. 



2 THE SUZAKU & SWIFT OBSERVATIONS 
OF MR 2251-178 

Suzaku (Mitsuda et al. 2007) observed MR 2251-178 in the 
XIS nominal pointing position for 287 ks between the 7— 10*'^ 
of May 2009. Data are included from the X-ray Imaging 
Spectrometer (XIS; Koyama et al. 2007) and from the PIN 
instrument of the Hard X-ray Detector (HXD; Takahashi 
2007), both of which were processed using version 2.3.12.25 
of the Suzaku data reduction pipeline. 

A complementary hard X-ray dataset obtained with 
the Swift Burst Alert Telescope (BAT), which observed 
MR 2251-178 as part of the 58 month all sky survey (Baum- 
gartner et al. 2010), is also included in our analysis. A sum- 
mary of the observations is included in Table [1] 
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Table 1. Summary of observation parameters 



Instrument 



Date 



Exposure 
(ks) 



Count rate 
(s-l) 



Flux 

(xlO^^^ ergcm"^ s~-^) 



Obs. ID 



Suzaku 



XlS-Fl 
XlS-Bl 
HXD/PIN 
HXD/GSO 



07/05/09 



136924 

103800 
89228 



2.037 ±0.003 
2.534 ± 0.004 
0.135 ± 0.003 
< 0.3 



4.29" 
4.35" 
5.48'' 
< 4.80'' 



704055010 



Swift 



BAT 



(14.5 ± 0.3) X 10~ 



6.63'' 



N/A 



" 2-10 keV flux. 

15-50 keV flux. 
" 50-100 keV flux 90% upper limit 

20-100 keV flux. 



2.1 XIS Data Reduction 

XIS data were selected in the 3x3 and 5x5 modes using 
ASCA grades 0,2, 3, 4 and 6. Standard XIS screening crite- 
ria were applied such that data were excluded if taken: (1) 
within 436 s of passage through the South Atlantic Anomaly 
(SAA), (2) within an Earth elevation angle (ELV) < 5°, 
and/or (3) with Earth day-time elevation angles < 20°. Hot 
and flickering pixels were removed from the XIS images us- 
ing the CLEANSIS script. Source spectra were extracted from 
within circular regions of radius 2.8', and background spec- 
tra were extracted from offset annuli of the same radius with 
care taken to avoid the corners containing the Fe^^ calibra- 
tion sources. The redistribution matrix flle and ancillary re- 
sponse flle were generated using the tasks xisrmfgen and 
XISSIMARFGEN, respectively. 

After checking for consistency (see i]3.1|l . spectra ob- 
tained from the two front-illuminated XIS and XIS 3 detec- 
tors were combined into a single source spectrum (hereafter 
referred to as XIS-FI) using mathpha in order to maximise 
signal-to-noise. Data from the back-illuminated XIS 1 (here- 
after XIS-BI) were not combined and were instead analysed 
separately. Only XIS data over the energy range 0.6-9.0 keV 
were included due to a degradation in S/N above 9.0 keV. 
Data were also ignored between 1.6-2.1 keV due to uncer- 
tainties associated with the the SiK edges intrinsic to the 
detector assembly of the XIS instrument. 

The background subtracted count rates during the ob- 
servation were 2.037±0.003 s"^ per CCD for the XIS-FI and 
2.534±0.004 s~^ in the case of the XIS-BI. The net exposure 
for the observation was 136.9 ks. All XIS spectra and corre- 
sponding response files were binned to sample the half-width 
half-maximum (HWHM) energy resolution of the detectors 
(i.e. '^GOeV resolution at 6keV). Counts were additionally 
grouped with GRPPHA to achieve a minimum of 50 counts 
per energy bin to enable the use of minimisation, which 
was used for all subsequent spectral fitting. 



2.2 HXD/PIN & GSO Data Reduction 

The HXD/PIN spectrum was extracted from the cleaned 
events files and was also processed according to the screening 
criteria described previously. The Non X-ray Background 
(NXB) was generated using the tuned event files made 
available by the HXD instrument team (Fukazawa et al. 
2009) with a count rate of 0.461 ± 0.001 counts s"\ The 



Cosmic X-ray Background (CXB) was simulated using the 
form of Boldt (1987) and resulted in a CXB count rate of 
~ 0.024 counts s^^ . Simultaneous good time intervals (GTIs) 
were found for both the source spectrum and the NXB using 
the MGTIME task and the NXB exposure time was increased 
by a factor of ten to reduce the effects of photon noise. The 
NXB and CXB were combined with equal weight using the 
MATHPHA operation and subsequently subtracted from the 
source spectrum within xspec. 

The background subtracted count rate in the HXD /PIN 
was 0.135 ± 0.003 counts s~^, corresponding to a 15-50 keV 
flux of 5.48 X lO"^'^ ergcm"^ s~^. The net exposure time was 
103.8 ks after detector dead time was taken into account 
with HXDDTCOR. Spectra were binned to a 5a level above 
the background per bin from 15-50 keV. 

For consistency we also reduced the data obtained with 
the HXD/GSO. The total (background subtracted) GSO ex- 
posure was 89.2 ks for a count rate of < 0.3 counts s"'^, which 
corresponds to 90% upper limit on the 50-100 keV flux of 
< 4.8 X 10"^^ ergcms"\ Note that the GSO data are not 
used during spectral fitting and are used solely to check the 
consistency of the time-averaged Swift/BAT dataset. 



3 BROAD-BAND X-RAY SPECTRAL 
ANALYSIS 

We performed a wide-band spectral analysis of MR 2251-178 
in the 0.6-180 keV energy range using xspec v. 12.6.0(/ (Ar- 
naud 1996) and heasoft v. 6.10. AU fits were modified by 
a Galactic column of AfH(Gai) = 2.42 x 10^° cm"^ (Dickey 
& Lockman 1990) which was not allowed to vary. We used 
a photoelectric absorber with Wisconsin cross-sections (the 
wabs component in XSPEC; Morrison & McCammon 1983) 
to model the Galactic absorption. Spectral parameters are 
quoted in the rest-frame of the source {z = 0.064) unless 
otherwise stated and errors are quoted to the 90% level 
for 1 parameter of interest (i.e. Ax^ = 2.71). The cross- 
normalisations for the XIS-FI and XIS-BI were allowed to 
vary in all fits and generally found to be within ±2 — 3% 
of each other. The HXD/PIN normalisation was set to 1.16 
of the XIS. The normalisation on the Swift/BAT was left 
as a free parameter throughout, and was found to be 4- 
5% greater than that of the XIS instruments. There was no 
strong variability (i.e. < 10% level) observed over the course 
of the observation. Lightcurves for the three XIS detectors 
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Figure 1. Lightcurves for the XISO (green triangles), XISl (red 
circles) and ttie XIS3 (black closed squares). Only small variations 
in the lightcurve were observed during the observations. Icr error 
bars are roughly the size of the plot points. 



are shown in Fig[T] In all spectral plots XIS data are binned 
to sample the half- width half-maximum energy resolution of 
the detectors (i.e. 60 eV at 6keV) unless otherwise stated. 



3.1 Initial Continuum Parameterisation 

We first considered the X-ray spectrum of MR 2251-178 in 
the 3-5 keV band of the XIS data. An initial fit of the 
continuum with a simple power-law of photon index F = 
1.56 ± 0.01 modified solely by Galactic absorption yielded a 
reasonable fit (with x^/degrees of freedom=212. 5/178). In- 
dividually, the photon indices for the three XIS spectra are 
all in good agreement with this value, and have indices of 
r = 1.55 ± 0.01, r = 1.55 ± 0.02 and V = 1.57 ± 0.02 for the 
XIS 0, XIS 1 and XIS 3 data, respectively. 

Extending the data to include the full 0.6-180.0 keV 
energy range reveals significant deviations from the sim- 
ple power-law fit and the fitting statistic is extremely poor 
(x^/dof = 4733.3/348). The residuals of this fit are shown 
in Fig. [2] (top panel). In the XIS data a clear soft excess 
can be seen at energies below around 0.7 keV and there is a 
noticeable discrepancy between 0.8-1.6 keV, presumably due 
to the presence of the warm absorber. The ubiquitous FeKa 
line is present at ~ 6.4 keV and the spectrum also appears 
to roll over above 10 keV in accordance with the presence of 
the high-energy cut-off reported by Orr et al. (2001). Note 
that as the Swift/BAT spectrum is time-averaged over 58 
months, we checked the HXD/GSO spectrum to determine 
the validity of the high-energy cut-off present in the BAT 
data. We find a 90% upper limit for the HXD/GSO flux of 
< 4.8 X 10"" ergcm"^ s"^ in the 50-100 keV band, which 
is consistent with the Swift/BAT flux over the same energy 
range (2.7 x 10~^^ ergcm~^s~^) and suggests that the roll- 
over is valid for the Suzaku observation. 

To parameterise the 0.6-180 keV continuum we initially 
flt the data with a model of the form wabs x (powerlaw -|- 
diskbb + zgauss) x zwabs, where the diskbb component is 
a simple parameterisation of the soft excess as an accretion 
disk blackbody, zgauss is a Gaussian used to model the 
Fe Ka fluorescence emission line and zwabs is a simple pho- 
toelectric absorber modelled in the source rest-frame (i.e. at 
z=0.064). 

With this model the fit statistic is drastically improved, 
but is still quite poor (x^/dof — 1111.1/340). The power- 
law component has a photon index of F = 1.67 ± 0.01 and 



is absorbed by A'^h = 1.1 x 10^ cm~ , which is similar to 
that found by Ramirez et al. (2008) for the lowly ionised 
absorber in the C/iandra/LETG observation of this source. 
The accretion disk blackbody used to parameterise the soft- 
excess has a temperature kT — 62^3 eV. The FeKa line is 
unresolved {a < 391 eV) and found at a rest-frame energy of 
6.43lg Q3 keV, which is consistent with an origin in neutral 
or mildly ionised material. As shown in Figure [2] (second 
panel), the poor flt statistic is due to the model being unable 
to adequately reproduce the spectral curvature between ~ 
2 — 6keV and in the hard X-ray band. This suggests that 
additional absorption is required. 

We therefore tested for a more complex absorption sce- 
nario by adding an additional neutral absorber with Ah = 



8.7: 



X 10 cm which covers 25 ± 1% of the source flux 



(modelled with zpcf abs in XSPEC). This gives a signiflcant 
improvement to the fit (x^/dof — 565.5/338 for two addi- 
tional free parameters) and the residual spectral curvature 
in the XIS data is no longer present. With a softer photon 
index of F = 1.83 ± 0.01, the new model was also able to 
replicate the observed fiux in the HXD/PIN data. Neverthe- 
less, visual inspection of the third panel in Figure [2] shows 
that the model is still unable to reproduce the curvature 
present above ~ 50keV (Figure [21 third panel). 

In order to flt the observed rollover in the Swift/BAT 
data we replaced the power-law with a cutof fpl component 
which models a power-law with an exponential high-energy 
cut-off (-Bcut)- The addition of the cut-off further improves 
the flt by Ax^ = 63.2 for 1 more free parameter and there are 
no longer any residuals in the BAT data (Figure O bottom 
panel). In this model F = 1.72±0.03 and Scut = HGtH keV, 
zwabs has Ah ~ (1.1 ± 0.1) x 10^^ cm~^, zpcf abs covers 
17 ± 2% of the source fiux with Ah = 9.9ti| x 10^^ cm"^ 
and the soft excess is parameterised by a kT = 57 ± 2 eV 
accretion disk blackbody. The final flt statistic is x^/dof = 
502.3/337 and is our best-flt continuum parameterisation. 

3.2 Absorption Lines 

Several absorption features are present below 3keV and in 
the FeK band. We modelled these absorption lines with 
Gaussian proflles flt in the rest-frame of the AGN. The 
statistical significances for each absorption line in this sec- 
tion were determined against the best-flt continuum model 
discussed above. Monte Carlo simulations (MC) were also 
performed to further assess the statistical signiflcances; the 
details of which are discussed in i|3.3l Table [2] shows a sum- 
mary of all line parameters. Detailed discussion regarding 
all line identiflcations, including the consideration of alter- 
native identiflcations, is presented in appendix A. 

3.2.1 The Soft X-ray band 

Below 2 keV two absorption lines are required (Fig. (3] top 
panel). The first is unresolved and detected at a rest-frame 
energy E = 0.77 ± 0.01 keV. This is consistent with the 
Fe i-xvi M-shell unresolved transition array (UTA) which is 
expected between 0.729-0.775 keV (Behar et al. 2001, Netzer 
et al. 2004). The line has an equivalent width (EW) of -7± 
2eV and its addition to the model improves the global flt 
by Ax^ ~ 65.6. The line is found to be >99.9% signiflcant 
from MC simulations. 
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Figure 2. Broadband data/model residuals of tlie Suzaku (XIS-FI: blaclc, XIS-BI: red, PIN: green open triangles) and Swift (BAT: 
blue open squares) data when fit with a variety of parameterising models ( ^3.11 1. From top to bottom: (1) Residuals when fit solely 
by the 3-5 keV best-fit power-law (2) Residuals of the same data when fit an accretion disk black body and power-law absorbed by a 
fully-covering neutral absorber. (3) As above, except with an additional partial-covering absorber included in the model. (4) As above, 
but with a cut-off power-law. Note that the Fe Ka has not been fit for illustrative purposes. 
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The second line has a rest energy of 1.29 ±0.01 keV and 
an EW=-4t5 eV. The line is highly significant (> 99.9% 
from MC simulations) and its addition improves the fit by 
Ax^ = 66.1. Unlike the UTA, this second line does not 
have an energy that corresponds to any expected strong 
X-ray transitions. There are several possible identifications 
for this line (see Appendix A) but the most conservative 
association which is internally consistent with our xstar 
model (35) is with the Fexxiv 2s— >-3p doublet which is 
expected at an mean energy of E — 1.165 keV. The line 
identification requires a blueshift of Vout = 0.11 ± 0.01c 
(33000 ± 3000 kms~^) assuming this identification. 

Two further Gaussian absorption lines are required at 
E = 2.52 ± 0.02 keV and E = 2.79 ± 0.03 keV in the rest- 
frame, with equivalent widths of —6 ± 2eV and —5 ± 2eV, 
respectively (Fig. [S] middle panel). Both lines are signifi- 
cant at the > 99.9% level from MC simulations and their 
addition improves the model by a further Ax^ = 15.9 
and Ax^ ~ 10.9, respectively. Theoretically, the strongest 
spectral features expected in the 2-3 keV energy range are 
those associated with the Is— ^2p transitions of S xv {E = 
2.461 keV) and Sxvi (E = 2.623 keV), but the rest-frame 
energies of the detected lines do not correspond to either 
of these which suggesting that they may be blueshifted. If 
identified with the S xv and S xvi transitions the measured 
line centroids indicate that they are blueshifted by Vout = 
0.02 ± 0.01c (6000 ± 3000 km s~^) and t^out = 0.07 ± 0.01c 
(21000 ± 3000 kms"^), respectively. Both lines are unre- 
solved and their widths were fixed to cr = 10 eV throughout. 



3.2.2 The FeK region 

There are further residuals in the data/model ratio present 
in the FeK region. In addition to the FeKa emission line 
there is a clear broad absorption feature in the 6.5-8 keV 
range, as can be seen in the bottom panel of Fig. O Note 
that this absorption is found in excess of the neutral Fe K 
absorption edge (expected at E=7.1keV) which is included 
as part of both the zwabs and zpcf abs absorbers. This ab- 
sorption trough was initially modelled with a broad Gaus- 
sian profile and its centroid energy left as a free parameter. 
The centroid energy and intrinsic width of the broad line 
are E = 7.57toll keV and a = 209twt ^V, and it has an 
EW of -26ti2 eV. The addition of the line results in a Ax^ 
improvement of 16.6 for 3 additional free parameters, and it 
is statistically required at the 99.3% level from MC simula- 
tions. With the addition of the broad absorption component 
the width of the Fe Ka emission line is now more tightly con- 
strained than before at o" < 117eV (while EW=25tg^ eV). 

In this energy range the spectrum is expected to be 
dominated by atomic features attributable to iron. In par- 
ticular, the Is— >-2p transitions are expected to be particu- 
larly strong, as are the many resonance transitions associ- 
ated with the Fe K edge near 7.1 keV (Kallman et al. 2004). If 
due to the strong Is— >2p transitions this immediately places 
a constraint on the possible identifications for the absorp- 
tion as being due to at least Fexvill {E ~ 6.5 keV), where 
the Fe ions become sufficiently ionised to have an L-shell 
vacancy for a Is electron transition. Above this ionisation 
state, the Is— >2p transitions occur between E ~ 6.5 — 6.6keV 
for Fexvill-xxiv, and at S ~ 6.7keV and E ~ 6.97 keV for 
FexxvHea and FexxviLya, respectively. If the broad pro- 





Fe XXV/XXVI 
(+ Fe K edge) 



Observed Energy (keV) 

Figure 3. XIS-FI (black open circles) and XIS-BI (red closed 
circles) residuals when fit with the best-fit continuum parame- 
terisation, as discussed in i|3.1l The positions of the absorption 
lines are clearly labelled. The XIS-BI data have been omitted in 
the bottom panel due to reduced S/N at FeK and the remaining 
XIS-FI data have been additionally binned by a factor of 3 above 
the HWHM energy resolution of the detectors for clarity. 



file is identified with Fe xxv Hea it requires a velocity shift 
of Vout = 0-iSto°o2 c (= 39000tS kms^^), while if identi- 
fied with Fexxvi Lya the velocity shift is slightly lower at 
Vout = 0.101^2 c (30000l^g°^ kms-i). 

An alternative interpretation is that the broad absorp- 
tion represents a blend of the above lines rather than a sin- 
gle discrete profile. To investigate the possibility that the 
profile is a blend of He- and H- like Fe we replaced the 
broad profile with two narrow (cr — 50 eV) absorption lines 
at fixed rest frame energies of 6.7 keV and 6.97 keV and left 
their common velocity shift as a free parameter. This fit 
was a slight improvement to that of the broad Gaussian ob- 
tained previously, and gives a Ax^ = 18.1 for 3 parameters 
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6 6.5 7 

Observed Energy (keV) 

Figure 4. Plots showing the FeK band. Top: the absorption 
feature fit with a broad Gaussian. Middle: the same absorption 
feature when fit as two separate Gaussian lines using the method 
outlined in i|3.2.2l 



of interest. Even so Fig. |4] shows that this model slightly 
under- predicts the absorbed flux of the Fexxv Hea line at 
E = 7.361^1 keV (rest-frame; E = 6.92 observed), while 
the Fexxvi Lya line is well fit. The resulting simultaneous 
velocity shift is v = 0.10 ± 0.01c with respect to the rest 
frame of the AGN. 



3.3 Montecarlo Simulations 

The Montecarlo simulations we used are analogous to those 
carried out by Tombesi et al. (2010a). Briefly: 

(i) We simulated S = 1000 XIS spectra between 0.6- 
10.0 keV based upon the baseline continuum parameterisa- 
tion model (outlined in i]3.1|l , using the fakeit command in 
XSPEC. Both the HXD/PIN and Swift /BAT spectra were 
not included in the simulation as they had no effect on the 
absorption line parameters. 

(ii) In each of the simulated spectra we stepped an in- 
verted Gaussian of a set flxed width (see below) between 
two pre-deflned energy bands in equal steps. After each step, 
the model was fit for the line normalisation of the Gaussian 
and the resulting overall Ax^ was recorded. To sample the 
entire XIS spectrum we performed this process twice, once 
to cover the FeK band and once to cover the soft X-rays. 
For the Fe K absorption we stepped a broad absorption line 
of a = 200 eV width every 100 eV between 4-9 keV, while 
in the soft X-ray band we stepped a narrow absorption line 
of 10 eV width every 25 eV between 0.6 and 4keV. Due to 
the XIS-BI having poorer S/N around FeK we used only 
the XIS-FI data in this band. Both the XIS-FI and XIS-BI 



were used in the soft X-rays where the S /N of the XIS-BI is 
higher. 

(iii) After S spectra were generated the associated grid 
of values then corresponded to how statistically likely it 
would be for a randomly occurring feature to have a Ax'^ 
improvement greater than a particular value. To assess the 
probabilities of the absorption lines detected, the Ax^ value 
we obtained from our Gaussian fits was then compared to 
this generated grid of values. If A'' of the simulated spectra 
had a feature of significance greater than this Ax^ value, 
the resulting detection confidence level of the measured line 
was 1 - N/S. 

The statistical probabilities for each line from MC sim- 
ulations are recorded in Table [2] 



4 SELF-CONSISTENT MODELLING 

With both the general continuum and the absorption lines 
parameterised, we proceeded to replace the simple parame- 
terising models (i.e. zpcf abs and zwabs) with a series of self- 
consistent model grids generated using the xstar photoion- 
isation code. Each xstar grid contains a series of models, 
each pertaining to a particular photoionised spectrum with 
a characteristic ionisation parameter {£,),& column density 
(Nu), and a redshift (z). In addition, each grid is attributed 
an intrinsic turbulent velocity (Kurb) upon generation which 
determines the intrinsic width of any emission and absorp- 
tion lines. The grid is then used as a simple multiplicative 
model in xspec to simultaneously fit the broadband contin- 
uum and any spectral lines which may be present. 



4.1 Fully covering models 

We first modelled the spectrum assuming that the pho- 
toionised absorbers fully cover the view to the ionising 
source, and that the soft excess represents an intrinsic com- 
ponent of the source continuum. We sequentially added XS- 
TAR grids to a baseline continuum model of the form wabs x 
(cutof f pl+ diskbb+zgauss) x xstar until an acceptable fit 
had been achieved. Solar abundances of Grevesse & Sauval 
(1998) are assumed throughout unless otherwise stated. 



4.1.1 Lowly ionised absorbers 

In the soft X-ray band two low turbulence (uturb ~ 
200kms~^) ionised absorption components are required to 
provide a good description to both the Fe UTA absorption 
line and the general spectral curvature below lOkeV (see 
Fig. El). The first absorber (Zone 1; Ax^ = 93.7 for 3 addi- 
tional free parameters) is lowly ionised and mainly respon- 
sible for the FeUTA, while he second absorber (Zone 2) is 
less significant to the fit {Ax^ = 36.9 for 2 more free param- 
eters) and mainly fits the curvature. These absorbers also 
contribute to a shallow Fe K edge in the Fe K band, however 
its depth is insufficient to adequately replicate the feature 
(see ^4.1.4|l . For fitting purposes the soft X-ray absorbers 
were assumed to be co-spatial and their common outflow ve- 
locity was fixed to Wout = 300 km s"'^ to be consistent with 
the UV outflows reported in this source by Monier et al. 
(2001). Allowing the outflow velocities to vary did not yield 
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Table 2. Summary of X-ray absorption and emission lines 



Line ID 


Brest (keV) 


Ei^h (keV) 


Flux" 


a(eV) 


EW (eV) 






MC 


FeUTA 


0.77 ±0.01 


0.729 - 0.775 


— 941 0^ 


10* 


-7± 2 


65.6 


> 


99.9% 


Fe XXIV 
(2s->3p) 


1.29 ±0.01 


1.165 


-38.olit;I 


< 44 




66.1 


> 


99.9% 


Sxv 


2.52 ±0.02 


2.461 


-17.2 ± 5.7 


10* 


-6 ±2 


15.9 


> 


99.9% 


Sxvi 


2.79 ± 0.03 


2.623 


iJ^-0_4.8 


10* 


-5 ±2 


10.9 


> 


99.9% 


FeKa 


"■^'-'-0.03 


6.4 


+14.0+«'^ 


< 117 


+25i^2 


52.9 






Fe xxv-xxvi 


7.5710-2 


6.7 and/or 6.97 


-8.9ir7 




99+15 
^^-11 


16.6 


99.3% 



* Denotes parameter was fixed at listed value 

° Line flux quoted in units of xlO~^ ergcm"^ . 

' Change in fitting statistic when adding a Gaussian to best-fit continuum model. 




Observed Energy (keV) 

Figure 5. Data/model residuals of the XIS-FI data when fit- 
ted with two lowly ionised soft X-ray absorbers. The position of 
the Fe L, S XV, S xvi and Fe K absorption is clearly apparent. For 
clarity, the XIS-BI data have again been omitted. Data between 
1.6-2.1 keV have been removed due to the instrumental Si K edge. 
See text for further details. 



any appreciable change in the fit statistic and they were thus 
left fixed throughout. 



4-1.2 Detection of a FeL-shell and Fexxv absorber? 

To model the absorption line at ~ 1.29 keV we used a higher 
turbulence (uturb ~ 5000 km s~^) grid to account for the ap- 
parent line broadening. Zone 3 (Ax^ = 74.0; see TablsO 
for parameters) fits both the 1.29 keV line (Fig|6} and also 
appears to fit the higher energy component of the FeK ab- 
sorption (Fig. at i? ~ 7.7 keV (rest-frame; 7.24 keV ob- 
served) . The outflow velocity found for this absorber is con- 
sistent with an identification of the FeK absorption with 
Fexxv Hea, however, the velocity is substantially larger 
than that found if the 1.29 keV line is identified as solely 
Fexxiv. This could be the result of the line at 1.29 keV be- 
ing a blend of Fexxill (2s^3p; expected at 1.127 keV) and 
Fexxiv (2s— !>3p). If this is the case, the outflow velocity 
would then be consistent with that found for the xstar grid. 
We hereafter denote this absorption component the "FeL- 
shell absorber" as at an ionisation of log^/ergcms~^ = 3.02, 
the most prominent ions are from Fexxill-xxv. 



4-1.3 Sulphur absorbers 

There is a level of ambiguity regarding the correct physical 
interpretation of the lines at 2.52 and 2.79 keV. In particu- 
lar, we were unable to obtain a satisfactory fit to the lines 
assuming solar abundances, and the lack of other distinct 
absorption lines meant that when fit with variable abun- 
dances, the required Sulphur over- abundance was entirely 
unconstrained. With the current data we therefore parame- 
terise the absorption lines using a column density comprising 
solely of Sulphur (A'^s). 

As discussed in H3.2.1I the most conservative interpre- 
tation is that the lines are identified with the Is— >-2p tran- 
sitions of He- and H-like Sulphur, respectively. In this inter- 
pretation the lines are well modelled by the addition of two 
further ionised absorbers. The first (Zone 4), which models 
Sxv at ~ 2.5 keV, is described log^/ergcms"^ = ^.AAILq^^ 
and a Sulphur column density of A^s = 1.91° ^ x lO^^cm"^. 
For a solar Sulphur/Hydrogen abundance ratio of 2.14x 10~^ 
(Grevesse & Suvaul 1998) this corresponds to a Hydrogen 
column of A^h = 8.9l'l3 * x lO^'^ cm~'^, which is similar to the 
column densities reported previously for this source (i.e. Orr 
et al. 2001; Kaspi et al. 2004; Gibson et al. 2005, Ramirez 
et al. 2008). The net outflow velocity of this absorber is 
Vout = 63001^200 kms-i (0.0211^;;°^;^ c), which is consis- 
tent with that found when fitting simple Gaussians. At this 
column density and ionisation parameter the grid also con- 
tributes He/3 (Is— >3p) and weak He7 (Is— >-4p) lines in ad- 
dition to strong Hea, however the energy separation of the 
lines is insufficient to fit the line at ~ 2.8 keV simultane- 
ously and a second grid is required to model the remaining 
residuals. 

The second Sulphur absorber (Zone 5; Ax^ = 32.9) 
which models the S xvi Lya is described by A'^s > 2.2 x 
10^'' cm"^ (= ATh > 1 X 10^^ cm"^) and log^/ergcms~^ = 
S.SltlH - The outflow velocity is found to be Vout = 21000± 
3000 kms~^ (0.07 ± 0.01c), which is again consistent with 
that found from simple Gaussian fltting. The fit statistic 
for the model assuming two zones of Sulphur absorption is 
xVdof = 352.2/316. 



4-1-4 ^ further highly ionised zone? 

While the soft X-ray and FeL absorbers contribute a shal- 
low FeK edge at ~ 7.1 keV (see middle panel of Figure 
[Tjl its depth is insufficient to fit the remaining absorp- 
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Figure 6. Ratio plots of the 0.8-1.6 kcV energy band after the 
inclusion of absorption zones 1 and 2, showing: Top: the resid- 
ual FeL line profile. Bottom: the same line but when zone 3 of 
the photoionisation model is added. The residuals are clearly no 
longer present. 
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Figure 7. Ratio plots of the Fe K band in the fully-covering XSTAR 
model. Top: when no ionised absorbers are included two residual 
profiles are apparent between 6.5-8.0 keV. Middle: when Zone 
3 is added the higher energy residual is no longer present. Bot- 
tom: when a further highly-ionised zone is added there are no 
remaining residuals. See text for further details. 



tion feature at £ ~ 7.36 keV (rest-frame, observed frame 
E ~ 6.92 keV). This feature can be tentatively (Ax^ = 8.0 
for 3 additional parameters; Fig [7] bottom panel) fit by the 
addition of a highly ionised and high column density ab- 
sorber (uturb ~ 5000 km s~^) which is outflowing with Vout = 
18000 ± 3000 kms"^ (0.06 ± 0.01 c). The parameters imply 
that the feature may be identified with Fexxvi (Is— >-2p), 
which is expected to be the must abundant ion at such an 



ionisation. This feature may be a higher velocity analogue 
of the outflowing Fexxvi absorption («out ~ 12700 kms~^) 
which was reported by Gibson et al. (2005) in a earlier Chan- 
dra/HETG observation of this source. 

^.1.5 Reflection Component 

The above absorbers in addition to the cutoffpl, diskbb 
and zgauss components provide a good description of the 
general broad-band continuum and the fit is marginally 
acceptable (xVdof = 352.2/316; nuU probability=0.08). 
However, for consistency we replaced the redshifted Gaus- 
sian modelling the FeKa emission with the ionised reflec- 
tion model reflionx (Ross & Fabian 2005), which self- 
consistently takes into account any reprocessing due to re- 
flection in a face-on system. Fixed Solar abundances were 
assumed and the photon-index of the reflected power-law 
was tied to that of the intrinsic. 

The addition of a lowly ionised reflection component 
(5 < 27ergcms^^) with a low reflection fraction [R < 0.2, 
where R = Q./2n and is the solid angle subtended by 
the reflector) reproduces the FeKa profile and improves 
the fit in the hard X-ray band, for a final fit statistic of 
xVdof = 342.79/316 (nuU prob=0.14), F = and 
Ecvit ~ 10lj;24 keV). The overall best-fit parameters for the 
fully-covering model are listed in Table [3] 

4.2 Partial covering models 

We also investigated the possibility that only a fraction 
/ < 1 of the line-of-sight source flux may be covered by the 
ionised absorbers. This would then require F for the power- 
law continuum to steepen to compensate for the increased 
absorption, and similar partial covering scenarios have been 
invoked to self-consistently model hard X-ray excesses in 
several AGN (e.g. NGC1365, Risaliti et al. 2009; 1H0419- 
577, Turner et al. 2009; PDS 456, Reeves et al. 2009). Mod- 
elling the hard X-ray excess in this way typically requires 
the intrinsic flux below lOkeV to be strongly absorbed by 
Compton thick material (i.e. A'h > 10^^ cm"'^). However, if 
the partial coverer is Compton thin, it may be able replicate 
the roll-over and soft excess. 

We thus constructed a low turbulence (uturb = 
200kms~^ partial covering model using the zxipcf model 
in XSPEC. Our starting model had the form wabs x 
(zxipcf xpowerlaw+reflionx) x xstar, where we have re- 
moved the diskbb component which previously parame- 
terised the soft-excess and reverted to a simple power-law 
with no high-energy cut-off. The overall best-fit parameters 
for our partial covering model are listed in Table 2] Im- 
portantly, this model still requires the presence of several 
fully covering absorbers to account for the FeL/FeK and 
Sulphur absorption, and these absorbers have best-fit pa- 
rameters which are consistent with those obtained in the 
fully covering model. A reflionx refiection component was 
also included for self consistency. In this geometry the par- 
tially covering clouds would need to be very compact in size 
and located close to the ionising source to only obscure a 
fraction of the emergent continuum. One possibility is that 
the partial coverer may be associated with a clumpy disk- 
wind where the covering fraction relates to the filling factor 
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[h] 



Table 3. Best- fit parameters for the fully covering model 



Continuum parameters 



Flux (erg cm ^ s ^) 
Luminosity (ergs~^) 



-Po. 5-2.0 
^2. 0-10.0 

io. 5-2.0 
.^^2. 0-10.0 



3.0 X 10'" 
4.5 X 10-" 
2.83 X 10** 
4.41 X lO** 



Best-fit model components 

Component Parameter Value 



Galactic absorption 
cutofEpl 

reflection 



diskbb 

Soft X-ray absorber 1 
(Zone 1) 



Soft X-ray absorber 2 
(Zone 2) 



FeL absorber 
(Zone 3) 



S XV absorber 
(Zone 3) 



Sxvi absorber 
(Zone 5) 



r 

Abund 



R 

rpc 
in 

norm 

"h 
logr 



2.42 X 10^" 
1.68 ±0.01 
10lli7 
1* 

< 27 

< 0.2 
62^5 

2.96+\ °i X 10^ 

(5.4 ±0.4) X 10^0 
-0 23+0'i0 



5.6 



/ 

"out 



logC-^ 

^out 



log?'^ 

^out 
"turb 



300 

200 

^I'l X 1021 

9 9-1+0.08 
^•^-^-0.07 

300 
200 

5.2tli X 1021 

o 09+0-08 
'-'■'-'^-0.03 

42000+3°°° 
5000 
8.9I37 X 10^1 

9 44+0.24 
^■**-0.22 



6300 



5100 
-4200 

1000 



log?" 

^out 
turb 



> 1.0 X 10^2 



3.37 



.+0.12 
19800+1°° 



1000 



77.2' 



114.6 



93.7 



36.9 



74.0 



76.4 



32.9 



" Absorber column density in units of cm'^. 
^ Cut-olf energy in units of keV. 

Accretion disk black body temperature, in units of eV 
' ionisation parameter in units of erg cm s"!. 
^ Outflow Velocity of absorber given in units of kms" 1. 
9 Assumed absorber turbulent velocity in units of kms'l. 

Signiflcance of component with respect to the final best-fitting 
model. 

* Statistical significance of REFLIONX component is given in 
comparison to a fit where Fe Ka is left unmodelled. 



0.1 



0.05 





redshift (z) 



-0.05 



-0.1 



Figure 8. Contour plot showing the redshift for the high- velocity 
soft X-ray absorber in the partial-covering model. 



of the denser clouds, while the more ionised absorbers ap- 
pear to fully cover the source along the line-of-sight. Such 
scenarios have been effective in explaining the broad-band 
spectral properties of objects such as MCG -6-30-15 (Miller 
et al. 2008) and 3C445 (Reeves et al. 2010), and are con- 
sistent with theoretical expectations that disk-winds are in- 
homogeneous outflows and contain clumps of dense, high 
column-density, material (e.g. Proga & Kallman 2004; Sim 
et al. 2010) 

In this case the partial covering model requires a very 
lowly ionised, high column density absorber which covers ap- 
proximately 30% of the source in order to fit the soft-excess 
and curvature below lOkeV (Table Q. As expected, the 
power-law is also much steeper (F — 1.96 ± 0.01) than that 
usually inferred for this source (F = 1.6 — 1.7). With just the 
partial coverer included significant residuals remained in the 
soft X-ray band and two low turbulence (uturb = 200 km s'^) 
fully- covering absorbers are also required to achieve an ad- 
equate fit (Zones 1 and 2 in Table |4)|. 

Interestingly, the outflow velocity of Zone 2 is Vout = 
27000 ± 3000 kms"^ (0.09 ± 0.01c), which is much faster 
than the 100-1000 km s~^ expected for typical warm ab- 
sorbers (e.g. Blustin et al. 2005; Crenshaw, Kraemer & 
George 2003). Nevertheless, the very high outflow velocity 
is strongly required in this model (see Figure (Sjl, and flxing 
the outflow velocity to the source rest-frame (i.e. Uout = 0) 
resulted in a worse fit by Ax^ = 46.6 (for one less free pa- 
rameter). As it is not required in the fully-covering scenario, 
the high outflow velocity of this component is likely to be 
highly model dependent and further observations at a higher 
resolution are required to determine its validity. 



^.2.1 The FeK band 

As in the fully-covering model, the reflionx component 
gives an excellent description of the Fe Ka emission line and 
the reflector is once again found to be both lowly ionised 
and have a low reflection fraction. The partial coverer pro- 
vides a very good fit to the residual absorption in the Fe K 
band (see inset of Figure. [9] bottom panel) and, in addition 
to the Fexxv Hea absorption which is being modelled by 
the FeL absorber, the lowly ionised partial coverer provides 
an excellent fit to the Fe K edge. 
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4-2.2 The high- energy roll-over 

While the overall fit statistic for the partial covering model 
is good (x^ = 335.6/308; null prob=0.13) the fit still 
appears to overestimate the flux at the highest energies 
(see Figure |9l bottom panel). We thus tried to replicate 
the roll-over by adding an exponential roll-over using the 
highecut model; however this only marginally improved the 
fit (Ax^ = 4.9 for one additional free parameter). Neverthe- 
less, the final fit statistic for the partial covering model is 
X^/dof = 330.7/307 (null hypothesis=0.168) which is com- 
parable to that obtained with the fully-covering model. 



5 DISCUSSION 
5.1 The soft excess 

We have demonstrated that the soft excess in MR 2251-178 
can be adequately modelled either as an intrinsic compo- 
nent of the emergent X-ray continuum (i.e. an accretion disk 
black body), or as the side effect of an absorption dominated 
(partially covered) spectrum below lOkeV. However, there 
are several further alternative interpretations which we also 
considered. 

5.1.1 Comptonisation 

A particular caveat of the fully-covering model is the man- 
ner in which we have parameterised the soft excess as be- 
ing the thermal emissions from a geometrically thin ac- 
cretion disk (diskbb), and thus as an intrinsic compo- 
nent of the X-ray spectrum. While simple and generally 
well fitting in a wide range of objects, such an inter- 
pretation typically requires a thermal continuum with a 
temperature higher than that expected for geometrically 
thin accretion disks (e.g. Shakura & Sunyaev 1974; Por- 
quet et al. 2004; Gierliriski & Done 2004). The expected 
temperature for a geometrically thin disk is approximately 
T(r) « 6.3 X 10'5(Macc/Medd)'/''M-'''*(r/i?,)-3/4K, where 
Ms is the accretion disk mass in units of 10* Mq . With 
(Macc/Modd) ~ 0.15 (see later) and Ms ~ 2.4 (Dunn et 
al. 2008) appropriate for MR 2251-178, and assuming most 
of the emission originates at the innermost stable circular 
orbit (risco = 3Rs for a Schwarzschild black- hole), we have 
kT ~ 20 eV. This is roughly 3 times cooler than the best- 
fitting diskbb temperature {kT ~ 62 eV) and implies that 
the soft excess cannot easily be attributed to thermal emis- 
sion from a standard accretion disk without modification. 

An alternative interpretation of the soft excess as an in- 
trinsic component is that it is the Wien tail of soft photons 
from the accretion disk which have been Comptonised by a 
hot corona located above the plane of the disk (e.g. Czerny 

6 Elvis 1987; Haardt & Maraschi 1991, 1993; Nayakshin, 
Kazanas, & Kallman 2000; Nayakshin 2001). If the distribu- 
tion of electrons in such a Comptonising corona is "quasi- 
Maxwellian" (i.e., has both a thermal and non-thermal com- 
ponent; see Coppi 1999), it may be able to produce both the 
hard X-ray power-law continuum for which it is typically in- 
voked, as well as the soft excess. 

To see whether this interpretation was feasible, 
we replaced the diskbb component in our best-fitting 
fully-covering xstar model with the comptt component 



Table 4. Best-fit parameters for the partial-covering model 



Continuum paramet^ 


ers 


Flux (erg cm-2s-M ^0.5-2.0 
^^2. 0-10.0 

Luminosity (ergs"-"-) ^0.5-2.0 
-^2. 0-10.0 


3.9 X 10"" 
5.3 X 10-" 
3.88 X 10** 
5.18 X lO** 


Best-fit model components 

Component Parameter Value 



Galactic absorption 
highecut 
power-law 
Reflection 



zxipcf 



Soft X-ray absorber 1 
(Zone 1) 



Soft X-ray absorber 2 
(Zone 2) 



FeL absorber 
(Zone 3) 



S XV absorber 
(Zone 4) 



S XVI absorber 
(Zone 5) 



Ecut 

r 

Abund 
R 

-^H , p GOV 

log ^pcov 
/gov 
^turb 

N- 

logr 



V' 



f 

out 
9 

turb 



logr 



V' 



f 

out 
9 

turb 



2.42 X lO^" 

> 244 
1.96 ±0.01 
1* 

< 25 

< 0.2 

(1.2 ±0.2) X 10^3 

r, r,7+0.02 
'-'•^'-0.03 

200 

(7.5 ±0.4) X 10^1 
1.79 ±0.02 

< 300 
200 

(2.llO;i) X 1021 
1.83 ±0.04 
27000 ± 3000 
200kms-i 



4.0 



,+0.5 
-0.6 



X 10^ 



loge 
^out 



loge 

/ 

out 
9 

turb 



6.4 



o r,i+0.23 
'^■"'^-0.03 

39000+*^°°° 
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^ Outflow velocity of the absorber given in units of kms^^. 

3 Assumed absorber turbulent velocity in units of kms^^. 

Significance of component with respect to the final best-fitting 
model. 

* Statistical significance of REFLIONX component is given in 
comparison to a fit where Fe Ka is left unmodellcd. 
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Figure 9. Unfolded best-fit models for (a) the fully- and (b) the partially-covering models. Inset in each panel is an enlarged view of the 
FeK region of each model. The la upper limit to the HXD/GSO flux is represented by the open circle in both plots and is consistent 
with the data obtained from the Swift/HAT 58month all-sky survey. Note that we have not included the tentative Fe xxvi component 
from the inset panel of (a). 



(Titarchuk 1994). The thermal temperature of the input 
accretion disk photons was assumed to kT ~ 20 eV, while 
the redshift was fixed to that of the source. The resulting 
fit was comparable (Ax^/dof = 335.1/308), and the other 
model parameters were all consistent with those obtained 
previously when using diskbb. The Comptonising electron 
temperature is kT < 5.3 keV and the optical depth of the 
corona is r ^ 0.40lo;3g. 



5.1.2 Ionised reflection 



Another alternative explanation for the soft excess in gen- 
eral is as a byproduct of ionised reflection from optically 
thick matter in the nuclear regions of the AGN, e.g. the ac- 
cretion disk, the broad line region (BLR) or the molecular 
the torus. However, due to the lack of a strong ionised reflec- 
tion component and the absence of a broadened iron line, 
this seems unlikely in the case of MR 2251-178. 
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5.1.3 Section summary 

In summary, the soft excess in MR 2251-178 is well modelled 
by either a thermal accretion disk black body or as the result 
of the Comptonisation of soft thermal photons from the ac- 
cretion disk. A partial covering interpretation of the dataset 
is also capable of replicating the soft-excess albeit with a 
steeper power-law photon index, i.e. V ~ 2.0, which is simi- 
lar to the value found commonly for radio-quiet quasars (e.g. 
Porquet et al. 2004). The lack of a strong ionised reflection 
component above lOkeV and the absence of a broadened 
FeKa line suggests than an origin through blurred reflec- 
tion is unlikely. Regardless of the true nature of the soft 
excess in MR 2251-178, the parameters that describe the 
fully-covering absorption components are largely invariant 
between the differing interpretations. 



5.2 The FeKa line 

From past observations of MR 2251-178 there has been sig- 
nificant uncertainty in the reported FeKa emission EW: 
i.e. 125^113^ eV {Gmga; Mineo & Stewart 1993), igolgf eV 
{ASCA; Reynolds 1997), 62^25 eV, {BeppoSAX; Orr et al. 
2001) 53 ± 20 eV {XMM-Newton; Kaspi et al. 2004) and 
23li3eV {Chandra; Gibson et al. 2005). 

A possible explanation for the uncertainty is that the 
previous observations lacked the spectral resolution and S /N 
in the FeK band in order to resolve the absorption feature 
blue-ward of the Fe Ka, and this may have artificially broad- 
ened the measured EW of the FeK emission line. Inciden- 
tally, we note that the two observations thus far which have 
reported the presence of Fexxv/xxvi absorption (this work 
and Gibson et al. 2005) have also reported the tightest con- 
straint on the Fe Ka EW. 

Even so, the FeKa EW we find for MR 2251-178 is stiU 
below average. For example, based on an archival study of 
36 sources in the Chandra/BETG archive (which includes 
MR 2251-178) Shu et al. (2010) have reported the narrow 
core of the Fe Ka emission line has a weighted mean EW 
of 70 ± 4 eV, which is substantially broader than that which 
we find here. Furthermore, from their large Suzaku study of 
nearby Seyfert galaxies, Fukazawa et al. (2010) have found 
that a mean Fe Ka EW of around 55 eV for sources which 
are absorbed by column densities similar to MR 2251-178 
(i.e. iVH<10^^cm-^). Given that MR2251-178 is a highly 
X-ray luminous AGN, the apparent weakness of the FeKa 
line could be a result to the 'X-ray Baldwin Effect' - an 
apparent anti-correlation between source X-ray luminosity 
and the Fe Ka EW, which has been well documented in the 
literature (e.g. Iwasawa & Taniguchi 1993; Page et al. 2004; 
Jiang et al. 2006; see also Shu et al. 2010 and Fukazawa et 
al. 2010). Furthermore, given that MR 2251-178 has a FRI 
radio morphology, it could also be due to a dilution of the 
observed X-ray spectrum by a component of X-ray contin- 
uum emitted in a radio-jet. The radio-loudness of an object 
is typically given in terms of the radio-loudness parameter, 
_Rl, defined by 



-Rl = log 



F 



(1) 



objects where _Rl > 1 are regarded as radio-loud, while those 
where i?L < 1 are radio-quiet. For MR 2251-178 Reeves & 
Turner (2000) found that i?L = —0.43 which suggests that 
any present radio-jet is very weak and unlikely to contribute 
in the X-ray band. 

A further possibility is that the weakness of the FeKa 
emission may be due to the inner accretion disk being 
shielded by the outfiow. If this is the case, this could also 
explain the lack of a strong ionised refiection component as 
the observed refiection would have to originate at distances 
further from the central engine, i.e. the BLR or the torus. 
Sim et al. (2010) have also recently shown that FeKa emis- 
sion can originate via reflection from dense material in an 
outflowing wind (with EWs~ 10s of eV), which may also 
be applicable here. Furthermore, the FeKa EW can also 
become suppressed in the presence of a sufficiently high- 
column absorber (e.g. Miller et al. 2010; Yaqoob et al. 2010) 
such as those found here. 



5.3 The multi-component outflow of MR 2251-178 

We have shown that the X-ray spectrum of MR 2251-178 can 
be adequately modelled with both a fully-covering and an 
absorption dominated partially-covered interpretation. We 
find that both of these models require the presence of a num- 
ber of fully-covering absorbers to account for the observed 
Fe UTA, FeL (Fexxiii-xxiv), Sxv, Sxvi and Fexxv-xxvi 
absorption lines. In this section, we turn our attention to 
these absorbers and investigate their geometries, kinematics, 
and, ultimately, estimate their energetic output. For simplic- 
ity we consider here the absorbers as per the fully-covering 
model, and focus on the Fe L absorber as it is likely to be the 
most energetic given its high Uout- For completeness, how- 
ever, the methods discussed here were individually applied 
to all absorbers in the fully-covering model and the results 
are summarised in Table [5] 

Note that as the parameters of the FeL absorber are 
largely invariant between the two continuum models, the 
results obtained here are applicable to either interpretation 
and are thus not model dependent. 



5.3.1 Distance 

The observed column density of a spherically symmetric ab- 
sorber a distance R from the ionising source is equal to the 
absorber density multiplied by the thickness of the shell of 
material; i.e. A^h ~ nSR. From the definition of the ioni- 
sation parameter, — Lion/nR^, where Lion is the ionising 
luminosity between 1-1000 Rydbergs (13.6 eV to 13.6 keV) 
and n is the hydrogen gas density, this then implies that the 
inner face of the absorber, R, is located at 



5R 



R=^ 



(2) 



where ohz is the monochromatic radio flux at 5 Ghz and 
F ° is the optical flux at A = 4400A. In this regime, 

4400A ^ ' 



from the central ionising continuum. Assuming that the ab- 
sorber forms a thin shell where 5R/R < 1, equation ((2)| then 
sets an upper limit on the distance to the inner face of the 
absorbing material, i.e. R < Lion/C-^H- From the best-fit 
continuum parameters of the fully-covering model. Lion is 
estimated to be of the order of 2 x 10*^ ergs"'^. Therefore, 
taking A^h ~ 4.4 x 10^^ cm^^ and log^/ergcms~^ ~ 3.3 
characteristic of the FeL absorber, we find R < 100 pc. 
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Table 5. Summary of absorption zone parameters 





Soft X-ray 1 
(Zone 1) 


Soft X-ray 2 
(Zone 2) 


Sxv 
(Zone 3) 


S XVI 
(Zone 3) 


FeL/FeK 
(Zone 5) 


Vout (kms 1) 

Reac (pc) 

Mout (Mgyr-i) 
b 

-Eout (ergs-1) 

-Eout(corr) (ergS^^) 


300* 
2 X 10^ 

20 
330006 

10-5 

~ 1045 6 
-10^0 


300* 
I uu 
20 
1006 
10-3 
~ 10*2 6 
~1039 


gOQQ+5100 

ouu 
0.06 
15006 
10-2 
~ 10'*«6 
~ 10*4 


198001-°° 

0.005 

6006 

10-3 
~ 10«6 
- 10« 


42000+3°°° 
1 on 
0.001 
27006 
10-3 
~ 10*86 
~ 10« 



'Denotes the parameter was fixed at this value. 



Assuming that the outflow escapes the system we can 
set a lower limit on the radial distance by considering the 
escape radius of the material. For a spherical geometry the 
radius at which material will be able to escape the BH is 
given by 



Ra. 



> 



2GM 



2c2 



(3) 



with _Rg = GM/iP' being the gravitational radius. Given 
that the FeL absorber has an outflow velocity of ~ 0.14 c, 
^csc 100 -Rg. This suggests that the outflow may have an 
origin consistent with an accretion disk wind or the BLR. 
Indeed, Kaspi et al. (2004) set a lower limit on the size of the 
BLR in MR 2251-178 of 0.04 ± 0.01 pc (f» 3000 Rg), which 
has significant overlap with the upper limit obtained here. 



5.3.2 Mass outflow rates 

We can also estimate the mass outflow rate. Assuming a 
quasi-spherical system: 



Mout ~ 4TvbR UenipV, 



4nb 



mpVout 



(4) 
(5) 



where 6 = Q/An < 1 is the geometrical factor allowing for an 
outflowing wind subtending sold angle f2. The parameters of 
the FeL absorber suggest that it is capable of expelling ~ 
27006 M0yr-^. For comparison, the accretion rate required 
to maintain the observed radiative bolometric output is: 



Ma, 



ibol 



(6) 



where Lboi is the bolometric luminosity and r) is the accre- 
tion efficiency. With i] = 0.06 for a Schwarzschild black- 
hole accreting at maximum efficiency, and adopting a UV- 
determined bolometric luminosity of 4.3 x 10^*5 ergs-^ (Dunn 
et al. 2008), we find that the AGN in MR 2251-178 requires 
an accretion rate of Mace ~ 1.3M0yr-^ to maintain its 
bolometric radiative output. This is orders of magnitude 
lower than the estimated mass outflow rate and implies that 
either, (1) the AGN will exhaust its fuel supply of accreting 
material over a relatively short timescale in order to main- 
tain an effective mass outflow rate much higher than Mace, 
or (2) the outflow is not spherical and has a small covering 



fraction or is rather clumpy (i.e. 6^1). Similar conclusions 
were also reached by Gibson et al. (2005) in their analysis 
of the Chandra/ WETS G observation. 

By assuming that Mout ^ Mace we can estimate the 
clumping/covering h factor of the outflow with respect to 
the central nucleus, i.e. 6 ~ Maec/Mout ~ 10-3. The result- 
ing inferred covering factors are listed in Table [5] We note 
that in principle 6 could be larger than this. However, this 
estimate allows us to set a conservative lower-limit to the 
energetics of the outflow in the next section. 



5.3.3 Energetics 

The kinetic luminosity of the FeL absorber is: 



1 . 2 

Eout ~ — Mout'^out^ ' 



10 erg s 



(7) 



which implies a conservative mechanical output of ~ 
10'*5 ergs-^ for h ~ 10-3 as discussed above. This makes the 
FeL absorber the most significant outflowing component in 
terms of overall energetics, despite its substantial putative 
clumping. The energetics of the other absorbers are listed in 
Tabled 

Assuming that the BH mass grows through accretion 
alone, we can then estimate a required 'active phase' for the 
QSO as the current estimated BH mass (2.4 X 10** Mq) di- 
vided by the accretion rate, i.e. ~ 10* yr. Combining this 
with the mechanical output calculated previously suggests 
that the FeL absorber has an estimated mechanical output 
of ~ 10®° erg over the lifetime of the AGN. This is compara- 
ble to the binding energy of a typical 10^^ Mq galaxy bulge 
and further supports the premise that such outflows may be 
energetically significant in AGN/Galaxy feedback scenarios. 



6 SUMMARY 

We have performed a broad-band spectral analysis of the 
radio-quiet quasar MR 2251-178 using recent observations 
from Suzaku and the Smft BAT. Below is a summary of our 
results: 

• The continuum shows considerable curvature above 
around ~ lOkeV and an apparent soft excess is present 
below IkeV. A weak {EW = 26lg'^eV), unresolved {a < 
117 eV), FeKa emission line from neutral material is present 
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at ~ 6.4 keV. The spectral curvature and soft excess can be 
equally well modelled as an artefact of an absorption dom- 
inated, partially covered, continuum or as a fully covered 
intrinsic continuum component. 

• Absorption lines due to the FeUTA and the FeL-shell 
are present below 2keV. The FeL line is consistent with 
being a blueshifted blend of 2s— >-3p transitions from Fo xxill 
and Fexxiv. Two further absorption lines are detected at 
~ 2.5 keV and ~ 2.8 keV, which we conservatively identify 
with Sxv Hea and Sxvi Lya. For this identification the 
lines require blueshifted velocities of ~ 0.02 c and ~ 0.07 c, 
respectively. From Montccarlo simulations we find that all 
soft X-ray absorption lines are significant at the > 99.9% 
level. 

• We also detect significant (99.3% from MC) absorption 
in the FeK band (E ~ 7.5 keV; rest-frame) which is consis- 
tent with the presence of blueshifted Fcxxv-xxvi. A single 
high turbulence {vturb = 5000 km s~^) ionised absorber pro- 
vides a good fit to both the FeL-shell line at ~ 1.29 keV 
and to a portion of the absorption present in the Fe K band, 
requiring an outfiow velocity of Vout ~ 0.14 c. 

• We derive distances to the absorbers and find that the 
Fe L absorption may originate at sub-parsec scales, possibly 
in an accretion disk wind. By balancing the mass outflow 
and mass accretion rates, wc estimate covering fractions (or 
'dumpiness') and subsequent energetics for all absorbers. 
The FeL absorber is consistent with being substantially 
clumped, by a factor of around 10~^, but still dominates 
in terms of the overall energetics with an estimated kinetic 
output of ~ 10*^ ergs"^ 
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APPENDIX A: LINE IDENTIFICATIONS 

The outflow velocities (and resultant energetics) of the ab- 
sorbers are contingent on the correct identifications of the 
absorption features. Most of the absorption features de- 
tected in the Suzaku spectrum have rest-frame energies 
which do not correspond to any strong features expected 
in the X-ray band, and have thus been conservatively iden- 
tified as being the blueshifted absorption from the strongest 
expected nearby spectral transition. Here we discuss alter- 
native identifications for these absorption lines. 



Al Fe XXIV 

The observed energy range of the suspected Fexxiv 
(2s^3p) line, i.e. 1.28-1.30 keV, is populated by an array of 
high-order transitions from low-abundance elements which 
makes a rest-frame identification unlikely. For example, the 
2s— >-nd transitions of Coxvill and Coxx are found in this 
energy range which can be ruled out based on very low as- 
trophysical abundances (i.e. Co/H ~ 10~*). Indeed, these 
transitions are also unlikely as they are expected to be very 
weak in comparison to their lower-order counterparts. 

The line at 1.29 keV is also unlikely to be identified with 
any blueshifted lines of the expected abundant elements (i.e. 
O, Ne) without an extreme outflow velocity. Indeed, if iden- 
tified with H- or He-like Neon (1.02 keV and 0.922 keV, re- 
spectively) the outflow velocity would be ~ 0.3 — 0.4 c. This 
is much larger than that typically found in AGN warm ab- 
sorbers (e.g. Blustin et al. 2005). Furthermore, an identifica- 
tion with the ls->2p lines ofNax/xi (1.12 keV and 1.24 keV) 
is unlikely as such Na absorption has not been reported in 
even the longest exposure AGN X-ray spectra (e.g. Kaspi et 
al. 2002). 

As an identification with any local Is— >2p transition is 
unlikely, we also considered nearby 2s— >-3p transitions. The 
most likely 2s— >-3p transition local to the observed ~1.29 keV 
line is that of Fexxiv at ~ 1.165 keV. For this identification 
the line would then require an outflow velocity of ~ 0.1 c. We 
note that while the oscillator strengths of the 2s— >3p lines of 
Fe XXV (expected at E 1.22 - 1.24 keV) are expected to be 
higher than those of Fe XXIV in this band, the required out- 
flow velocity assuming this identiflcation {vout 0.04 — 0.06 c) 
would be incompatible with that obtained from xstar flt- 
ting (see i|4.1.2|) . 



A2 Sxv/xvi 

The lines at E = 2.52 ± 0.02 keV and E = 2.79 ± 0.03 keV 
do not correspond with any expected transitions in the rest- 
frame. The only transitions which have consistent rest-frame 
energies are those associated with extremely low abundance 
elements which are not expected to be observed in the X-ray 
spectrum (e.g. Kr, Rb, Mo). 

The higher energy line is consistent with the Is— i-Sp 
He-like Sulphur line at i5 ~ 2.88 keV but such an identifica- 
tion is again unlikely given that at the best-fitting ionisation 
stage (Table [S]) H-like Sulphur is expected to be the domi- 
nant ion (Kallman & McCray 1982). Furthermore, if the line 
was caused by the He-like Sulphur ls^-3p line, we would still 
expect to see much stronger Is— >2p transitions, which are 
not present. Thus a Vout = origin appears to be ruled out 
for both of the lines and they are most likely be blueshifted. 

The most likely ions in the 2-3 keV energy band which 
are expected to be strong in the X-rays are He- and H-like 
Sulphur at E = 2.461 keV and E = 2.623 keV, respectively. 
An identification with Silicon (e.g. Sixiv Lya at 2.01 keV) 
would require a substantially larger soft X-ray outflow ve- 
locity V ~ 0.2 — 0.4 c and thus seems unlikely; as does an 
identification with any lower Z ions for similar reasons. 
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A3 Fexxv/xxvi 

While their oscillator strengths suggest that the absorption 

features between 7-8 keV are most likely Fexxv/xxvi, sev- 
eral other elements are also expected to have transitions 
above ~ 7keV which may interfere with this identifica- 
tion. Most notable of these are H- and He-like Nickel which 
have Is— ^2p lines at rest-frame energies of ~ 7.8 keV and 
~ 8.1 keV, respectively. The observed energy of the absorp- 
tion feature is consistent with neither of these transitions 
without significant redshift, however, making an identifica- 
tion unlikely. An identification with lower ionisation Nickel 
could be plausible but the column densities at which mod- 
erately ionised Nickel would be observed are unacceptably 
high (i.e. A^h > 10^'' cm~^) and with an abundance ratio of 
Fe/Ni~ 20 (Grevesse & Suvaul 1998) we would also expect 
to observe hugely dominant Is— >^2p Fe transitions, which are 
not observed. Thus an identification with Ni would appear 
to be ruled out. 

For the tentative high-ionisation Fe K absorber an iden- 
tification with any ions other than Fe xxvi is unlikely given 
the high ionisation as Fexxvi is expected to be the only 
abundant ions at the best fit ionisation parameter (Kallman 
et al. 2004). Furthermore, while Zone 5 may be contribut- 
ing Fexxiii-xxiv Is— >-2p lines to the FeK band, Fexxv is 
expected to be the dominant ion at the measured ionisation 
parameter and we adopt it as our most probable line iden- 
tification for the contribution of Zone 5 to the FeK band. 
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